Entanglement generation between distant atoms by Lyapunov control 
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We show how to apply Lyapunov control design to the problem of entanglement creation between 
two atoms in distant cavities connected by optical fibers. The Lyapunov control design is optimal in 
the sense that the distance from the target state decreases monotonically and exponentially, and the 
concurrence increases accordingly. This method is far more robust than simple geometric schemes. 



I. INTRODUCTION 

Atoms, or their artificial counterparts, quantum dots, 
in cavities or traps have great potential for applica- 
tions in quantum communication, metrology and in- 
formation processing. Since entanglement is a crucial 
resource in quantum computation and communication, 
the preparation of maximally entangled states is a cru- 
cial task. Nonlocal interactions between two physi- 
cal qubits are required to generate entanglement and 
there have been numerous proposals to effect such in- 
teractions, especially for atoms trapped in distant cavi- 
ties 5, 6, 7, 8, 9, 10, 11, 12, 13^14], and similar 
schemes are conceivable for artificial atoms such as quan- 
tum dots. Some of these proposals make use of continu- 
ous feedback in open quantum systems pH] but most are 
based on Hamiltonian systems, and in most cases only 
simple geometric control schemes are employed to cre- 
ate the maximally entangled state. These methods have 
the advantage of simplicity but unfortunately often suffer 
from robustness issues. 

In this work we explore an alternative control design 
inspired by Lyapunov functions [T5 H T61 fTf l [181 [T9l [20 | 
[5TJ [2U [531 to find control designs for robust en- 
tanglement creation. Lyapunov control design has the 
advantage of being sufficient simple to be amenable to 
rigorous analysis, and much is known about their con- 
vergence properties, robustness and stability. In partic- 
ular such design can be shown to be highly effective for 
systems that satisfy certain sufficient conditions, which 
are roughly equivalent to the controllability of the lin- 
earized system [501 [25| . Unfortunately, this appears to 
be a strong requirement not satisfied by many physical 
systems. However, in certain cases, in particular for sys- 
tems like the two-atom model proposed by Mancini and 
Bose [121 , we can circumvent these restrictions by consid- 
ering the dynamics on a subspace and successfully apply 
Lyapunov control to create maximally entangled states 
from certain initial product states in robust fashion. 

The paper is organized as follows. In Sec.|TT]we briefly 
review the distant-atom model and the geometric con- 




FIG. 1; Two cavities Ci and C2, each of which contains a two- 
level atom, are connected into a closed loop through optical 
fibers. The off-resonant driving field A generates an effec- 
tive non-local Hamiltonian Hsf / while the two local resonant 
lasers generate the local Hamiltonian Hiocai- 

trol scheme proposed in |13j to generate entanglement. 
In Sec. jlllj we briefly review Lyapunov control and show 
how to apply it to the problem of steering the system 
from certain product state to one of the four Bell state in 
a robust fashion. We will consider two control paradigms: 
one is to control the local Hamiltonian which is easier 
to implement experimentally; the other is to control the 
non-local interaction Hamiltonian, which might be pos- 
sible for certain systems. 



II. TWO-DISTANT-ATOM MODEL AND 
GEOMETRIC CONTROL 

We consider a two-qubit model where the qubits are 
encoded in two atoms or two quantum dots in distant 
cavities connected into a closed loop by optical flbers, as 
illustrated in Fig.jl] It was shown in [13j that eliminating 
the interacting light fleld between the two atoms in the 
dispersive regime leads to an effective Hamiltonian for the 
two- atom system of form H^ot — ^^locai + ^^cff, where the 
local Hamiltonian induced by interaction with resonant 
light and the effective interaction Hamiltonian are: 

ifiocai = B{X ®I + I®X) (la) 
H^s^2JZ^Z (lb) 
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where X, Y , Z are Pauli operators and / is the iden- 
tity operator, and the coupling constant B = 77 J where 77 
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interaction time t (1/J) 



FIG. 2: Concurrence as a function of the interaction time for 
the geometric control scheme for different values of coupling 
B. Achieving unit concurrence requires requires switching off 
the Hamiltonian precisely, due to the fluctuations concurrence 
curve. If the control Hamiltonian is switched on or off too 
early or too late, even by a small amount, the concurrence of 
the final state may be reduced significantly. The three sub- 
figures on top of the main figure shows the zoom-in of the 
plots. 



should be sufficiently smaller than 1 to ensure the deriva- 
tion of Hcf[ remains valid. 

This Hamiltonian can be used to generate a maximally 
entangled state from the initial ground state by turning 
on TJtot for a critical time to before switching the field 
off [13] . Broadly speaking, by applying a constant Hamil- 
tonian we effectively perform a rotation about a fixed axis 
in the two-qubit space, and with the correct timing we 
can choose the rotation angle such as to ensure that the 
system state ends at the correct target state. However, 
plotting the concurrence of the final state versus the in- 
teraction time (Fig. [2]) shows that achieving very high 
fidelity with respect to the maximally entangled state re- 
quires very precise switching as the concurrence is subject 
to small fluctuations. In the model we have assumed a 
fixed coupling strength J and controllable local field B. 
We see that increasing B significantly reduces the time 
required to prepare a maximally entangled state but also 
increased the magnitude of the fluctuations. E.g., for 
B — 0.1 the fluctuations around the peak are only about 
1% but it takes 157 time units to reach a maximally en- 
tangled state. For B = 0.4 on the other hand, we can 
prepare a maximally entangled state in about 1 /8 of the 
time but the concurrence fluctuations increase by a factor 
of approximately 15. Therefore, although such design for 
entanglement generation is quite simple, it is not robust 
against imperfections of switching time. 



III. LYAPUNOV CONTROL DESIGN 

In the previous section we have seen that the method 
of entanglement generation by switching a constant field 
on for a fixed amount of time is highly sensitive to small 
switching time errors. Ideally, we would like a control 
scheme where the concurrence of the two qubits converges 
to 1 asymptotically, and better without any fluctuations. 
In that way, the control is robust against switching time 
errors. A simple method that seems well suited to this 
task is Lyapunov-based design. Roughly speaking, the 
idea of Lyapunov control is to choose a suitable so-called 
Lyapunov function V and then try to find a control that 
ensures that V is monotonically decreasing along any dy- 
namical evolution. 

In the time scale where the Hamiltonian evolution is 
still a good approximation, many physical systems satisfy 
the the quantum Liouville equation (with h = I) 

p=-i[Ho + f{t)Hi,p], 

where we have assumed the Hamiltonian has two parts: 
Hq is the system Hamiltonian and Hi is the interac- 
tion Hamiltonian, with the interaction coupling constant 
modulated by the function f{t). For example, for a two- 
level structure of a single atom, Hq = is the energy 
splitting, and Hi = is the dipole interaction between 
the laser and the atom, with a varying /(t) by modu- 
lating the laser amplitude. The fact that /(t) can be 
varied is very crucial from control point of view, since 
this degree of freedom allows us to design the dynamics 
to derive the desired evolution. 

We can define a general control task thus: for a given 
target state p^, for example, a maximally entangled state, 
we wish to find a control function /(i), such that the 
system state p(f) will converge to pd, as i — > oo. In many 
applications, we allow Pd(€) to evolve under iJo, and the 
control requirement becomes p(t) Pd(f) as t — *■ oo, 
which is generally known as tracking control |25j . In the 
following we assume: 

Pd = -i[i?o,/Od] 

Motivated from the theory of Lyapunov function and 
the Hilbert Schmidt distance ||jo(t) — Pdii)\\2, we define 

y{p.Pa) = \\\p- Paf = \Tr\{p- pdf]. (2) 
Assuming k > 0, if we choose 
]{t) = !{p(t),pd{i)) = K^x{pd{t)\-iHi,p(t)\), (3) 
we find that for Y{i) = V{p{t),pd{t)), 
V{t) = -f{t)i:r{pd{t)[-iHi,p{t)]) - -Kf{tf < 0. (4) 

Hence V is a Lypunov function and the value of V mono- 
tonically decreases along any solution {p{t) , pd{t)) . More- 
over, every solution {p{t), Pd{t)) converges to an invariant 



set E, called the LaSalle invariant set, on which V van- 
ishes. 

Discussions on Lyapunov-based design in terms of den- 
sity operators have been analyzed [HJ In par- 
ticular, many target states can be shown to be almost 
globally asymptotically stable if the Hamiltonian satis- 
fies certain demanding conditions: (i) Hq be strongly 
regular and (ii) Hi be fully connected ^31. The for- 
mer condition translates into the requirement that Hq 
have distinct transition frequencies between any pair of 
energy levels. This rules out systems with degenerate 
or equally spaced energy levels. The latter condition is 
even more demanding. In the basis where Hq is diagonal, 
all the off-diagonal elements of Hi must be non-zero, i.e. 
transitions between any two energy level of Hq can be 
realized. When the strict conditions on the Hamiltonian 
do not hold, for most cases, the target state can be shown 
to be no longer asymptotically stable, and we no longer 
have p{t) Pd(t), implying that the control design be- 
comes ineffective. This really restricts the applicability 
of the method especially for higher-dimensional systems, 
including two-qubit models and spin chains. 

However, for high-dimensional systems with Hamilto- 
nian not satisfying the above conditions, it is still possible 
to make the target state asymptotically stable on a sub- 
space, where the Lyapunov control can be applied effec- 
tively. In the following, for the two-distant-atom model 
(Fig. [T|, we illustrate how the Lyapunov control design 
can be utilized to drive the system state from a product 
state to a maximally entangled state, despite the fact 
that the full Hamiltonian of the system clearly does not 
satisfy the strict conditions set out above. 



IV. LYAPUNOV CONTROL DESIGN FOR 
ENTANGLEMENT CREATION 

For the two-distant-atom model with Hamiltonian ([T|), 
we can either choose the control Hamiltonian Hi to be 
the local Hamiltonians Hi = i/iocai or the effective cou- 
pling Hamiltonian Hi ~ HcS, depending on which sce- 
nario is easier to implement for a particular physical sys- 
tem. 



A. Local Control 

First, let us consider controlling the local Hamiltonian. 
In this case we choose Hq = H^s = 2J{Z ® Z) and Hi — 
H\ocai = vJ{^ (E) I + I <E) X). To make the Hamiltonian 
easier to analyze, we transform from the Z-eigenbasis 
{|0),|1)} to the X-eigenbasis {|+),|— )}■ In this basis, 
the matrices for the Hamiltonian are rewritten as 
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FIG. 3: Local Control: control fields obtained from Lyapunov 
design for different values of k and distance between the sys- 
tem state and the Bell state l^"*"). The control design is ro- 
bust in that the field amplitude gently decreases to zero, and 
the semilog distance plot shows that the convergence to the 
target state is not only monotonic but also exponential with 
the converging rate determined by k. 



and it is easy to see that the eigenvectors of Hq are the 
Bell states 
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To generate maximally entangled state, we can choose 
Pd = |<J>"'")($+|, for instance, and the control f{t) = 
KTr{pd[—iHi, p{t)]), according to ([3|. Notice that Hq 
and Hi do not satisfy the strict condition in Section |III[ 
Thus this design cannot drive every state to the target 
state, but we can see that if the initial state of the system 

is p{0) ^\ + +){++\oT \ ) ( I then the state will 

converge to the target state. In fact, in the Bell-state 
basis, the Hamiltonian can be written as 



Hq = 2 J 



/I 

10 

0-10 

VO -1/ 



Hi = 2r; J 



^0 ON 

10 

10 

vO Oy 



where 



Hq = 2 J 



/O IN 

10 

10 

Vl Oy 



Hi = 2t]J 



'10 





vO -ly 



pd = 



^0 ON 

10 



^0 Oy 



P(0) 







ON 

1 ±1 
±1 1 

Vo Oy 



4 



For states initially prepared in the subspace S spanned 

by I + +) and | ), we clearly see that the dynamics 

under the Hamiltonian Hq + f{t)Hi will be confined in 
that subspace, and thus we can consider the dynamics on 
this two-dimensional subspace iS where the Hamiltonians 
and state take the form: 



Hq = 2 J 



Hi = 277 J 
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The results in [53] now guarantee that all solutions in 
S except for |<i>~) will converge to the target state. The 
control field varies smoothly and steers the system gently 
to the target state as shown in Fig. [3] Moreover conver- 
gence is exponential and 



\f{t)\^K\i:r{i[p{t),pa{t)]Hi)\ 

^K\\l[p{t),pd{t)]Hl\\ 

<K\\i[p{t),pS)]\\-\\Hi\\ 



(6) 



shows that /(<) is bounded and we can choose k to ensure 
that |/(t)| is sufficiently small and the approximations 
inherent in the model remain valid. 

The method can also be utilized to increase the entan- 
glement in the initial state, i.e, to prepare a maximally 
entangled state starting with a partially entangled one. 
More specifically, if the system initial starts in the state 
IV'o) = All -I- -I- A2I ) then the control design pro- 
duces a control field that steers the system from this state 
to the desired maximally entangled state Choosing 
Pii = |<i>~)(<i>~| instead, we can similarly prepare |<i>^), 
and it can be verified that steering the state to l^*") sim- 
ply requires inverting the sign of the control field. Thus, 
not only can we prepare a maximally entangled state, but 
we can select which state we prepare. 

If the coupling constants of the local Hamiltonian for 
the two atoms are not exactly identical, e.g., if i?iocai — 
r}J{X iS) I + kl (E) X) then changing to the X-basis gives 
Hi = r/Jdiag(l + k,l — k, —1 + k, —1 — k), which trans- 
forms to 



Hl^TjJ 



/ l-k\ 

1+k 

1+k 





(7) 



Thus for k ^ 1 we can also steer the system from the 



product states 



or I to the Bell state |<i>*), i.e. 



for this two-atom model Lyapunov control can be used 
to prepare any of the four Bell states. 

One limitation of the scheme is that the initial state 
must be in the subspace S, for example, S = span{| -|- 

+), I )}, for the control to be effective. This is not 

a shortcoming of the proposed control scheme, however, 
because we can see from the structure of Hq and Hi 
that the control system is decomposable, hence not con- 
trollable on the whole space [26]. More specifically, the 
dynamics on the orthogonal subspaces S and S'^ are inde- 
pendent, and subspace populations are conserved quanti- 
ties. Thus, for the above Hamiltonian, no control exists 
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FIG. 4: Interaction control: control fields obtained from Lya- 
punov design for different values of k and distance between 
the system state and the target state with pd{0) = |^'''){^''^|. 
The control design is robust in that the field amplitude gently 
decreases to zero, and the semilog distance plot shows that 
the convergence to the target state is not only monotonic but 
actually exponential, although unlike in the local control case, 
\\p{t) — pd{t)\\ does not converge to 0. The final p{t) is still 
maximally entangled with unit concurrence. 



that steers population from subspace S to 5^ and vice 
versa. 



B. Interaction control 

Instead of controlling the atoms locally, we can alter- 
natively control the nonlocal Hamiltonian Hcs, if the un- 
derlying physical system allows. In this case we choose 
Hq = 'qJ{X®I + I®X) and Hi = 2J{Z®Z), or in the 
X-eigenbasis 



Hq = 277 J 
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The Bell states are no longer the eigenstates of Hq. 
Hence, for pd(0) = |<i>+)(<i>+|, the target state is also 
evolving with time, but for p{Q) = | + +) the dynamics 
is still confined to the subspace S spanned by | -I- and 

I ). Therefore, the dynamics can again be reduced to 

a 2D subspace on which we have 



Hq = 277 J 
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as well as 



Pd(0) = 



1 /I 1 



2 \1 1, 

where the orbit of pd{t) is the equator of the Bloch sphere. 
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FIG. 5: Evolution of concurrence C under Lyapunov control 
for different values of k for local control (top) and interaction 
control (bottom) shows monotonic convergence to 1. Insets 
show the error, i.e., 1 — C decreases effectively exponentially. 



From the analysis in |25] we can conclude that all so- 
lutions in S will converge to the equator of the Bloch 
sphere, i.e., states of the form 



1 



1, 



which corresponds to the LaSalle invariant set E of the 
original problem satisfying 



f I e-*"^ 





V" 1. 



Thus, we can no longer guarantee p{t) — > Pd{t) as 
t — !■ +00, i.e., that the state converges to a particular 
Bell state. This is illustrated in Fig. |4] which shows that 
the distance from the target state still decreases mono- 
tonically and exponentially but the asymptotic value of 
y(p(t), Pd{t)) for i — > oo now depends on k and is gener- 
ally larger than zero. However, since all the states in the 
set to which p{t) converge are maximally entangled, we 
can still steer the system to a maximally entangled state, 
and the concurrence still increases monotonically to one 
(see Fig. |5]) but the relative phase a of the state we con- 
verge to now depends on the exact initial state and the 
feedback strength k. 

Strictly speaking, as the control f{t) reduces to zero, 
the norm of Hiocai will cease to be significantly smaller 
than that of i?off, rendering the approximations made 



in the derivation of H^s invalid, unless we reduce the 
strength of i?iocai accordingly. However, in practice, 
the system should already have reached a state with 
significant entanglement before the model becomes 
invalid. 

V. CONCLUSION 

We have shown how to apply Lyapunov control to the 
problem of generating entanglement between two distant 
two-levels atoms in cavities connected by optical fibers. 
Given the Lyapunov control design, despite the fact that 
the sufficient condition for a target state to be asymptot- 
ically stable is not satisfied on the whole state space, we 
can still ensure it is almost globally asymptotically stable 
on certain subspace. Therefore, within that subspace we 
can drive the system from a product state to a maximally 
entangled state. The Lyapunov control design has the 
advantage of much greater robustness compared to sim- 
ple geometric schemes, and optimality in the sense that 
the distance from the maximally entangled target state is 
monotonically decreasing, and the convergence speed is 
exponential. We have discussed two control paradigms: 
to control the local Hamiltonian, as well as to control the 
effective interaction Hamiltonian between the two atoms. 
In both cases we can generate a maximally entangled 
state from an initial product state: for the formal case 
the system state will converge to a stationary Bell state, 
while for the latter case the relative phase of the final 
state will keep varying under the Hamiltonian, since the 
target state is non-stationary. Moreover, in the latter 
case, the model becomes invalid when the control am- 
plitude becomes sufficiently small. Therefore, the former 
control paradigm is preferable. The Lyapunov control de- 
sign can be also used to steer partially entangled states 
to a maximally entangled state, however, the control is 
only effective for initial states in the subspace where the 
target state is asymptotically stable. This is not a limita- 
tion of the control design, however, but a consequence of 
the fact that the controlled system is decomposable into 
two orthogonal subspaces on each of which the dynamics 
is invariant. In this sense, the Lyapunov control design 
is as effective as is possible within the constraints of the 
model. 
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